function y = computeDamageFunctionDerivative(element,etilde) %compute F'(etilde)
Dnn = element.Dnn;
g = computeTractionLaw(element,etilde);
dg = computeTractionLawDerivate(element,etilde);
y = Dnn/g - Dnn*etilde*dg/g^2;
end
